Ising transition driven by frustration in a 2D classical model with SU(2) symmetry 
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We study the thermal properties of the classical antiferromagnetic Heisenberg model with both 
nearest (Ji) and next-nearest (J2) exchange couplings on the square lattice by extensive Monte 
Carlo simulations. We show that, for J2/J1 > 1/2, thermal fluctuations give rise to an effective Z2 
symmetry leading to a finite-temperature phase transition. We provide strong numerical evidence 
that this transition is in the 2D Ising universality class, and that Tc — > with an infinite slope when 
J2/J1 1/2. 
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Since the milestone papers by Hohenberg and by Mer- 
min and Wagner, Q it is known that in two-dimensional 
systems a continuous symmetry cannot be broken at any 
finite temperature, and only systems with a discrete sym- 
metry can show a finite-temperature phase transition. In 
this regard, a proposal by Chandra, Coleman and Larkin 
(CCL) :2\ opened a new route to finite-temperature phase 
transitions in systems with a continuous spin-rotational 
invariance: CCL argued that the presence of frustrat- 
ing interactions can induce non-trivial discrete degrees of 
freedom, that may undergo a phase transition at low tem- 
peratures. In particular, in Ref.i2, the authors considered 
the antiferromagnetic Heisenberg model with both near- 
est (Ji) and next-nearest neighbor (J2) couplings: 

W = Ji ^ S, • Sj + J2 ^ s, • , (1) 

n.n. n.n.n. 

where Si are spin S operators on a periodic square lat- 
tice with N = L X L sites. For J2/J1 < 1/2, the classical 
ground state has antiferromagnetic Neel order with pitch 
vector Q = (ti", tt), while for J2/J1 > 1/2, the classical 
ground state consists of two independent sublattices with 
antiferromagnetic order. 3] The ground state energy does 
not depend on the relative orientations between the mag- 
netizations of the two sublattices, and the ground state 
has an SU(2) x SU(2) symmetry. Following Henley's anal- 
ysis of the XY case, CCL showed that both quantum 
and thermal fluctuations are expected to lift this degen- 
eracy by an order by disorder mechanism ^ and to select 
two collinear states which are the helical states with pitch 
vectors Q — (0, tt) and (tt, 0) respectively, reducing the 
symmetry to SU(2)xZ2. CCL further argued that the 
Z2 symmetry should give rise to an Ising phase transi- 
tion at finite temperature and provided an estimate of 
the transition temperature. 

The interest in this model raised recently with the 
discovery of two vanadates which can be considered as 



prototypes of the J1—J2 model in the collinear region: 
Li2VOSi04 and Li2VOGe04. I10 Indeed, although the 
value of J2/J1 is not exactly known, |3,.8ii.S>J all estimates 
indicate that J2 ^ Ji- Li particular, NMR and muon spin 
rotation magnetization in Li2VOSi04 provide clear evi- 
dence for the presence of a phase transition to a collinear 
order at Tc ~ 2.8K. While several additional ingredi- 
ents, like inter-layer coupling and lattice distortion, [loj 
are probably involved in the transition, the basic expla- 
nation relies on the presence of the Ising transition pre- 
dicted by CCL. 

However, CCL's predictions have been challenged by a 
number of numerical studies. Using Monte Carlo simu- 
lations, Loison and Simon reported the presence of 
two phase transitions for the XY version of the classi- 
cal model for J2/J1 > 0.5: A Kosterliz-Thouless transi- 
tion, as expected for XY models, followed by a transition 
which is continuous but does not seem to be in the Ising 
universality class since their estimates of the critical ex- 
ponents depend on the ratio J2/J1. More recently, the 
8=1/2 Heisenberg case has been investigated by Singh 
and collaborators [l^ using a combination of series ex- 
pansion methods and linear spin-wave theory. They show 
that, if there is a phase transition, it can only occur at 
temperatures much lower than that predicted by CCL for 
S=l/2, and they argue that Tc is actually equal to zero. 

In this Letter, we show, on the basis of extensive Monte 
Carlo simulations, that the classical limit of the model of 
Eq. Q , where spins are classical vectors of length 1 , in- 
deed undergoes a continuous phase transition at a finite 
temperature, that the critical exponents agree with the 
Ising universality class, and that, modulo minor adjust- 
ments of CCL's estimate, Tc is in good quantitative agree- 
ment with CCL's prediction in the range Ji < J2 where 
their approximation is expected to be valid. However, 
contrary to CCL's prediction, we show that Tc goes con- 
tinuously to zero when J2/J1 — ^ 1/2, and we argue that 



2 



this is due to a competition between Neel and coUinear 
order at finite temperature in this parameter range. 

Starting from the original spin variables Si, we con- 
struct the Ising-like variable of the dual lattice: 



(S, - Sfc) ■ (Sj- - S;) 

|(s.-Sfc).(s,-soi 



(2) 



where {i,j,k,l) are the corners with diagonal {i,k) and 
(j, I) of the plaquette centered at the site a of the dual 
lattice. The two collinear states with Q = (tt, 0) and 
Q = (0,7r) have CTq, = ±1. It is important to stress 
that the normalization term does not affect the critical 
properties of the model. It is only introduced to have 
a normalized variable. The Ising-like order parameter is 
defined as M„ = (l/N) J2a ^a- 

We have performed classical Monte Carlo simulations 
using both local and global algorithms as well as more re- 
cent broad histogram methods [I^ (details will be given 
elsewhere Q|) to calculate the temperature and size de- 
pendence of several quantities including the Binder cu- 
mulant, the susceptibility and the correlation length as- 
sociated to Ma, as well as the specific heat, for sizes up 
to 200 X 200 and for several values of J2/J1 between 1/2 
and 2. For reasons discussed below the critical behaviour 
is easier to detect for small values of J2/J1, and we first 
concentrate on J2/J1 = 0.55. 

As a first hint of a phase transition, we report the 
temperature dependence of the susceptibility defined by 
X = {N/T){{Ml) - (|M^|)2) for different sizes [see 
Fig. n^a)] . If there is a phase transition, this suscepti- 
bility is expected to diverge at Tc in the thermodynamic 
limit, and indeed, the development of a peak around 
T/Ji = 0.2 upon increasing the system size is clearly 
visible. To get a precise estimate of Tc, we have calcu- 
lated Binder's fourth cumulant of the order parameter 
defined by: Ui(T) = 1 - {W^) ji^Mlf . This cumulant 
should go to 2/3 below Tj, and to zero above Tc when the 
size increases, and the finite-size estimates are expected 
to cross around T^. Binder cumulants for different sizes 
are reported in Fig. db) , and they indeed cross around 
T/Ji ~ 0.197. In Fig.[Ii;c), we report U^'yT) as a func- 
tion of 1/L for several temperatures around 0.197. Ex- 
cluding temperatures for which [/4 clearly increases or 
decreases with 1 /L leads to the remarkably precise esti- 
mate Tc/Ji = 0.1970(2). 

To identify the universality class of the phase transi- 
tion, we have looked at the finite-size scaling of several 
quantities. The critical exponents v and 7 can be ex- 
tracted from the dependence of the peak position of the 
susceptibility Tc{L) = Tc + a x i^^/" and from its value 
x{L,Tc) ~ L''/'^ as a function of L. Using the value 
of Tc deduced from Binder's cumulant, the fits lead to 
1/ = 1.0(1) and ^jv — 1.76(2). A more precise estimate of 
the exponent v can be obtained from the temperature de- 
pendence of the second-moment correlation length ^ |l5l | , 
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FIG. 1: (a) Finite-size susceptibility x[L)/L and (b) Binder 
cumulant U4, as a function of the temperature for different 
sizes; (c) Binder cumulant U4, as a function oi\/L for different 
temperatures (0.1965,0.1966, . . . , 0.1973 from top to bottom). 
Lines are guides to the eye and the horizontal line marks the 
value for C/4 at the critical point for the Ising universality 
class. All data are for J2/J1 = 0.55. 



extracted from the Fourier components of the correlation 
functions of the Ising-like variable 0. By considering 
only values such that ^ < L/6, where the finite-size ef- 
fects are found to be negligible, it is possible to have a 
very accurate value of the critical exponent from the fact 
that, for T > T^, = A{T - T^Y . In Fig. 13 we report 
the behavior of the correlation length ^ as a function 
of the temperature. By performing a three-parameter 
fit (for A, Tc and v) we obtain T^/Ji = 0.1965(5) and 
f = 1.00(3). This value of Tc is compatible with the 
estimation given by Binder's cumulant. 

These exponents agree with those of the Ising univer- 
sality class in 2D (1/ = 1 and 7 = 7/4). A cross-check 
for this universality class comes from the measure of the 
critical exponent a, related to the divergence of the spe- 
cific heat per site, Cma.^{L) ~ L"/". Indeed, the value of 
a is the fingerprint for the 2D Ising universality class, for 
which we have a = and a logarithmic divergence of the 
specific heat as a function of L. In Fig. |2| we show the 
results for the specific heat per site. We have obtained a 
very accurate fit of the maximum of the specific heat per 
site as a function of L with the known expression for the 
leading finite-size corrections of the 2D Ising model, [l^ 
Cmax(T) = aQ + ai log(L) -I- a2/L, with ao, ai and 02 fit- 
ting parameters, while a power law is clearly inadequate. 
These results are consistent with a = 0. 

Finally, if the phase transition is indeed Ising, Binder's 
cumulant at Tc should reach the universal value U^iTc) ^ 
0.6107 in the thermodynamic limit. The non- 

montonous behaviour of Binder's cumulant with the size 
prevents a precise extrapolation, but this value is not 
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FIG. 2; Inverse of the correlation length as a function of 
the temperature for different sizes of the lattice and J2/J1 ~ 
0.55. The critical exponent u and Tc can be extracted from 
the behavior of as a function of the temperature. The 
arrow marks the resulting Tc- 
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FIG. 3: (a) Specific heat per site C as a function of the 
temperature for different sizes of the lattice and J2/J1 = 0.55. 
The lines are giudes to the eye. (b) Maximum of the specific 
heat per site Cmax(i) as a function of L. The line is a three- 
parameter fit (see text). 

incompatible with our numerical data [see Fig. ^c)]. 

Altogether, these results show unambiguously that 
a phase transition occurs for J2/J1 = 0.55 at Tc — 
0.1970(2) and give strong evidence in favour of 2D Ising 
universality class, in agreement with CCL's prediction. 
The same analysis can be repeated for different values 
of the frustrating ratio J2/J1, and the complete phase 
diagram is shown in Fig. ^ where we report as a 
function of J2I J\- While the critical behaviour is ev- 
erywhere consistent with Ising, it turns out that finite- 
size effects become more and more severe upon increasing 
J2/ Ji , preventing a meaningful determination of Tc with 
available cluster sizes beyond J2/J1 — 2. Indeed, for 
large ratios J2/J1, physical quantities such as the sus- 
ceptibility and the specific heat only exhibit broad peaks 
while the mean value of the order parameter goes very 
smoothly to zero, a behaviour typical of strong finite-size 
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FIG. 4: Monte Carlo results for the critical temperature Tc 
as a function of the frustrating ratio J2/J1. The line is an 
extrapolation of the large J2 data down to J2 = (see text). 

effects. This we believe can be traced back to the width 
of the domain walls between domains with Q — (tt, 0) 
and Q = (0, tt), which we have estimated by studying 
systems with fixed boundary conditions. Details will be 
presented elsewhere, [1^ but the width increases very 
fast with J2/J1, starting around 10 lattice spacings for 
J2/J1 ^ 1/2, and already reaching values of the order of 
40 lattice spacings for J2/J1 — 1-5. Since cluster sizes 
should be significantly larger than the width of the do- 
main walls to observe the critical behaviour, we are not 
able to go beyond J2/T1 ~ 2. 

Let us now discuss the dependence of Tc upon J2 / Ji 
(see Fig. 0J|. Two regimes can clearly be identified: (i) 
A large J2 regime (J2/J1 > 1) where Tc/Ji scales more 
or less linearly with J2 / Ji ; (ii) A smaller J2 regime de- 
fined by (J2/J1 — 1/2) ^ 1 where Tc vanishes with an 
infinite slope as J2/J1 1/2. For J2/J1 ~ 1/2, this 
disagrees with CCL's prediction that Tc/J2 reaches its 
maximum when J2/J1 — > 1/2. However, since their ap- 
proach is based on an expansion in J1/J2, this is not 
a final blow, and CCL's predictions should be tested in 
the large J2 regime, where the approximations are bet- 
ter controlled. CCL's central criterion for estimating Tc 
is the equation Tc ~ E(Tc){S,N(Tc)/a)^ , where EiTc) is 
the energy barrier to go from one domain to the other 
through the intermediate canted state where sublattice 
staggered magnetizations are perpendicular, and ^n{T) 
is the Neel correlation length of each sublattice. To get 
a quantitative estimate of Tc based on this approach, we 
have solved this equation using the exactly known tem- 
perature dependence of ^n{T) for the classical antifer- 
romagnet [iq and a corrected expression for E{T). 
This leads to Tc = 0.768 J2/(l + 0.135 ln(J2/Ji)). The 
best way to check this approach would be to detect the 
logarithmic correction, but unfortunately this would re- 
quire to go to temperatures much larger than what we 
can reach, and in the temperature range available to our 
simulations, CCL's prediction reduces to Tc ~ 0.77 J2. 
This prediction is in good agreement with our results: 
In the large J2 regime, Tc indeed scales linearly with 
J2, within error bars our high J2 data extrapolate to 
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at J2 = 0, and the slope is equal to 0.55. Note that the 
slight difference in slopes is not significant since including 
a constant factor in front of E{Tc) in the self-consistent 
equation for would change the slope. So altogether 
we believe that the present results support CCL's pre- 
diction when J2/J1 is not too close to 1/2. In that re- 
spect, we note that a similar estimate of Tc can be per- 
formed on the basis of CCL's criterion for S=l/2 using 
recent estimates of £,n{T) for the S=l/2 Heisenberg an- 
tiferromagnet on the square lattice, j2flj which leads to 
= 0.496 J2/(l -t- 0.78 ln(J2/Ji)). These estimates are 
considerably lower than those used in Ref. IT^ and with 
these new estimates we believe that a phase transition 
cannot be excluded on the basis of the numerical results 
of Ref. ^3 In fact, although the precise values of Ji and 
J2 in LiV0Si04 are still a matter of discussion, [7l|3,|9j we 
note that the order of magnitude of these new estimates 
is consistent with the experimental results of Ref. 0- 

Finally, let us discuss the behaviour of Tc close to 
J2/J1 = 1/2. The disagreement with CCL suggests 
that one ingredient is missing in that range. In fact, we 
found that the fluctuations have a clear Neel character 
when J2/ Ji ~ 1/2, implying that, although the coUinear 
phase is energetically favoured at zero temperature when 
J2/J1 > 1/2, thermal fluctuations favour the Neel state, 
as quantum fluctuations indeed do. We then expect a 
crossover to take place between a high-temperature Neel 
phase and a low temperature coUinear phase. To esti- 
mate the crossover temperature, we start from the low- 
temperature expansion of the free energy per site for clas- 
sical spins which reads, when only quadratic modes are 
thermodynamically relevant, 
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where eo is the ground-state energy per site and Uq are 
the frequencies of the quadratic modes. For the present 
purpose, the coefficients of this expansion should be de- 
termined in the limit J2/ Ji — > 1/2 from below and above 
for the Neel and coUinear states respectively, and the 
cross-over temperature Tq is determined by /nccI = f col- 
li turns out that, in the hmit J2/J1 — > 1/2, X^gl^w, 
has the same value for the Neel and coUinear phases, so 
that the linear term drops from this equation, leading 
to To oc {J2/J1- 1/2)1/2 gjjj^g gNcci ^ _2j^ ^ 2J2 and 
^0°'^ = — 2J2.I23I Now, the Ising transition can only oc- 
cur below To since the system should already have short- 
range coUinear fluctuations. Then, when the intrinsic 
temperature scale of the Ising transition as determined 
by CCL is larger than Tq, we expect the transition tem- 
perature to be of the order of Tq. This argument thus 
predicts that, as J2/J1 — > 1/2, Tc should vanish with 
an inflnite slope and with an exponent equal to l /2.|22j 
This is in qualitative agreement with our numerical data, 
which clearly indicate that the slope is infinite, and are 



consistent with an exponent in the range 0.3 — 0.5. Large 
scale numerical simulations are in progress to try to reflne 
this estimate. 

In conclusion, we have established the presence of a 
flnite-temperature phase transition in the classical J1 — J2 
Heisenberg model on a square lattice when J2/J1 > 1/2, 
we have provided strong arguments in favour of Ising 
universality class, and we have determined Tc as a func- 
tion of J2/J1 with high accuracy, showing in particular 
that it vanishes with an infinite slope when J2/J1 1/2. 
These results, together with revised estimates of Tc in the 
S=l/2 case following CCL, are expected to set the stage 
for further theoretical and experimental investigations. 
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